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Abstract 

We compute power-suppressed corrections to the 777 and r]'^ transition form factors Q'^F^[n')-i{Q'^) 
\ arising from the end-point regions x — > 0, 1 by employing the infrared-renormalon approach. The 

\ contribution to the form factors from the quark and gluon content of the r/, r/' mesons is taken 

00 ' into account using for the rj — rj' mixing the SUf{3) singlet r]i and octet 773 basis. The theoretical 

■ predictions obtained this way are compared with the corresponding CLEO data and restrictions 

O . on the input parameters (Gegenbauer coefficients) figC^i)) -^2(^1)' ^^ivs) in the distribution 

amplitudes for the rji, r/g states with one nonasymptotic term are deduced. Comparison is made 
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with the results from QCD perturbation theory. 
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I. INTRODUCTION 



The electromagnetic transition form factors (FFs) -^Ma(S Lpf lightpseudoscalar M = 
7r°, 1], 1]' mesons were the subject of much theoretical jH Q, M B M S B B EOl and 



experimental research in recent years. For instance, the CLEO collaboration reported 
about rather precise measurements of the 7r°7, 777 and r]''-/ transition FFs that stimulated 
interesting theoretical investigations aiming to account for the obtained experimental data 
within the framework of QCD. One of the key objectives of such analyses is to model 
the 7r°, T] and t]' mesons distribution amplitudes (DAs) and, in the r], r( case, to extract 
information on their gluon components. 

Indeed, it is known that the physical 77 and 77' mesons can be represented as superpositions 
of a flavor SUf{3) singlet rji and octet rjs state 

\r]) = cos9p\r]s) - sin6'p|?7i) , 

\t]') = sm9p\r]s) + cos6'p|?7i) . (1) 

Unlike the octet rjg state, the SUf{3) singlet r/i contains a two-gluon component fl2], which 
even absent at the normalization point /Xq , appears in the region > /Xq owing to quark- 
gluon mixing and renormalization-group evolution of the rji state DA. The rj and t]' mesons 
(cf. Eq. dD)) receive a gluon contribution due to the gluon content of the rji state. Because the 
meson-photon transition at leading order (LO) is a pure electromagnetic process, the gluon 
components of the 77 and rj' mesons can contribute directly to the 777 and 7/'7 transitions 
only at next-to-leading order (NLO) due to quark box diagrams. They also affect the LO 
result through evolution of the quark component of the f], rj' meson DAs. Contributions to 
the 7/7 and r]''y transition FFs, originating from the gluon content of the rj and 7/' mesons, 
were recently computed within the framework of the standard hard-scattering approach 
(HSA) of the perturbative QCD (PQCD) and estimates of the expansion parameters in the 
meson DAs were given. 

The gluon contributions to the 777 and 7/'7 electromagnetic transition FFs are subdom- 
inant. But in some exclusive processes, like the B meson two-body nonleptonic exclusive 
and semi-inclusive decays, which involve the rj and rj' mesons, their gluon contribution can 
potentially play an essential role in explaining the experimental data (see Ref. and refer- 
ences cited therein). The reason is that in these processes the gluon component of the rj and 
rj' mesons contributes to corresponding hard-scattering amplitudes already at LO of pertur- 
bative QCD. Hence, the gluonic parts of the rj, rj' meson DAs, deduced from the 777, ?7'7 
data, are important input ingredients in studying a wide range of exclusive processes, given 
that they are universal, i.e., process- and frame-independent quantities. 

The HSA and the perturbative QCD factorization theorems 3|, at asymptotically large 



values of the momentum transfer Q^, lead to reliable predictions for exclusive processes. But 
in the momentum-transfer regime of a few GeV^, experimentally accessible at present for 
most exclusive processes, power-suppressed corrections (l/Q^)", n = 1,2,3 .. . may play an 
important role in explaining the experimental data. In order to evaluate such corrections, 
the QCD running- coupling (RC) method, combined with the infrared (IR) renormalon ap- 
proach, was proposed [H, B 3, S El] • This method allows one to evaluate power-behaved 



contributions in exclusive processes arising from the end-point regions x — 0, 1. In this 
manner, power corrections to the electromagnetic FFs FuiQ"^) (M = tt, K) [3, 113 ; to 
the transition FFs Fu-yiQ"^) (M = 7r°, rj, rj') [l[ B], as well as to the gluon-gluon-7/' ver- 
tex function [l^ were computed. Power corrections can also be obtained by means of the 



2 



Landau-pole free expression for the QCD coupling constant [17[. This analytic approach 
was used to calculate in a unifying way power corrections to the electromagnetic pion FF 
and such to the inclusive cross section of the Drell-Yan process jl8. lol . 

Power corrections to the 777 and Vj^ electromagnetic transition FFs within the RC method 
were computed in Ref. and predictions for the structure of the DAs of the r] and 77' mesons 
were made. In the present work we extend this sort of investigation by also including into 
the calculation of the 777, ?7'7 transition FFs the power corrections originating from the 
gluonic content of the 77, 7/' mesons that were not taken into account in Ref. This will 
enable us to extract their DAs from comparing our theoretical predictions with the CLEO 
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data 

The paper is structured as follows: Sec. II contains the required information on the hard- 
scattering amplitudes for the 7717 and 7/87 transitions, accounting also for the gluon content 
of the r/i state. The DAs of the SUf{3) singlet rji and octet rjg states are considered and their 
evolution is taken into account. In Sec. Ill we compute the 777 and 1]''-/ transition FFs within 
the RC method and obtain the Borel resummed expressions for them. The asymptotic limit 

— >■ cxD of these FFs is explored and the standard HSA leading-twist predictions for the 
FFs are recovered. In Sec. IV we perform a numerical analysis and compare our results with 
the CLEO 11| data with the aim to extract constraints on the rj and 7/' meson DAs. Finally, 
Sec. V contains our concluding remarks. 



II. SUf (3) SINGLET AND OCTET COMPONENTS OF THE r/7, r]'-f TRANSITION 
FORM FACTORS 

The meson-photon electromagnetic transition FF FujiQ'^) can be defined in terms of the 
invariant amplitude F'^ of the process^ 

7*(gi)+7fe)^M(p) (2) 

in the following way 

F'^ = 7e2FM^(Q')e^^"^e,(g2)giag2/3 , (3) 

where ey{q2) is the polarization vector of the real photon and = —qf- The FFs of the 
777 and 77'7 transitions are sums of the corresponding singlet F^^{Q^) and octet F^^{Q^) 
contributions 

FM,m = F^,m + F^.iQ') ■ (4) 

The FF of the octet state F^^{Q^) and the quark- related component of the FF of the 
singlet state, Ffjj^((5^), can be computed by employing the results obtained for the pion- 
photon transition FF [i^]. The latter is known at 0{as) of PQCD [2^. More recently, 
also a part of O^a"^) corrections were computed ^|. The gluonic component of the singlet 



S / 

contribution was rust recently calculated within the framework of the HSA of the 

perturbative QCD in Ref. p. 

In accordance with the PQCD factorization theorems, at large momentum-transfer, the 
FFs F^^{Q^) and Fy^^i^Q"^) can be represented in the form of a convolution of the corre- 
sponding hard-scattering amplitudes with the quark and gluon components of the DAs of 



Hereafter M denotes the 77 or rj' meson. 
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the rji and r]s states, 



Qs(/iR) 

47r 



and 



(5) 



47r 



(6) 



where all quantities above are renormalized, i.e., are singularity- free, and the symbol 
denotes the convolution 



T^{x) ® (f){x) = [ dxTii{x)(f){x) 
Jo 



Here the functions T^ q{x) and ^^(x, Q^, /il) the hard-scattering amplitudes for the 
partonic subprocess 7 + 7* — » g + g at LO and NLO, respectively, and ^(x, Q^, is the 
NLO hard-scattering amplitude for 7 + 7* ^ g + g, with /i^ being the factorization and 
renormalization scales. In Eqs. (0) and (jH)), /m are the M- meson decay constants, Cp = 4/3 
is the color factor, and A^i and are numerical constants, each depending on the quark 
structure of the associated rji, rjs states 



el - lei 



(7) 



The hard-scattering amplitudes q(x), T^j ^(x, Q^, and i(x, Q^, are well- 
known and are given by the following expressions 
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(8) 



where x = 1 — x. 

The next ingredients needed for computing the FFs and are the meson- 

decay constants and the distribution amplitudes ^^^(x, 0^^(x, and /^|) 
for the r^i, T^g states. The decay constants are defined as the matrix elements of the 



axial- vector currents J^g with z = 1,8 



(9) 



In the octet-singlet basis the constants can be parameterized by two methods. One is to 
follow the pattern of state mixing (cf. Eq. [T])) 



fl = /s cos Op , = -fi sin Op , 
= fs sin 9p , f^, = fi cos 9p , 
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(10) 



where the decay constants /i, /g and 6p are given by 2^ 

h = i.nu, fs = i.26U, ep = -i5A° (11) 

with = 0.131 GeV being the pion decay constant. 

The second method employs a two-mixing-angles parametrization: 

fri = fs COS 6^8 , 4 = -/i sin Oi , 

/«=/8sin^8, /^i=/icos^i (12) 
with the mixing angles 9i and 6'8 provided by ji^] 

^1 = -9.2°, ^8 = -21.2°. (13) 

Parametrization (jlOp leads to simple expressions for the physical FFs in terms of F^^^(Q^) 
and F^8^(Q^); viz., 



F^^(Q2) = F^^^(Q2) cos^p - F^,-,(Q') sin a 



F^.^(g2) = F^3^(Q2) sin^p + F^i^(Q') cos^p , (14) 

where the form factors F^j^((5^) and F^g^((5^) are determined by expressions (0) and ©, 
but with the decay constants fl^ replaced by /j. In our numerical computations we shall 
use both schemes the conventional one-angle mixing scheme and also the two- mixing- angles 
parametrization. 

The main question still to be answered concerns the shape of the DAs of the r/i and rjg 
states. In general, a meson DA is a function containing all nonperturbative, long-distance 
effects, which cannot be calculated by employing perturbative QCD methods. Nonetheless, 
as a direct consequence of factorization, the evolution of the DA with the factorization scale 
/i| is governed by PQCD. Input information at the starting point of evolution, i.e., the 
dependence of the DA on the variable x at the normalization point /Xq, has to be extracted 
from experimental data or derived via nonperturbative methods, for example , Q CD sum 
rules with nonlocal condensates j2^ (see also 2^1) or instanton-based models 2^| at some 



(low) momentum scale, characteristic of the particular nonperturbative model. 

Because of mixing of the quark-antiquark component with the two-gluon part of the DA, 
the evolution equation for the DA of the flavor singlet pseudoscalar r]i state has a 2 x 2 
matrix form fl2\. The solution of this equation is given by the expressions 



/ 



[x, /ip) = 6xx 



1+ E {Bl 



and 



V 



6^(x, /ip) = XX 



n=2A.. 



n=2A.. 



«s(/io) 



' rn n 



«s(/io) 



X — X] 



(15) 



/So 



Qs(/^o) 



Cn-l(^ ^) 



(16) 



Here C^/^(z) and Cy^(z) are Gegenbauer polynomials. Detailed information concerning the 
parameters p^, and the anomalous dimensions 7", 7" can be found in Ref. [l^. In Eqs. 



5 



(fT^ and (Uni) the coefficients and will be considered as free input parameters, the 
values of which at the normalization point /Iq determine the shapes of the DAs 0''(x,/ip) 
and 

In our calculations we shall use a phenomenological DA for the rji state containing only 
the ffist Gegenbauer polynomials 6*2^^ (x — x) and Ci^'^{x — x) (i.e., i?! 7^ 0, -Bf 7^ and 
59 = 5^ = for all n>2) 

C2^'^{x — x)= 6(1 — 5xx), Ci^'^{x — x) = 5{x 
Under this assumption, the DAs assume the following forms [isf 



x] 



(17) 



()^^(x, /ip) 



Qxx 



1 + A(;Uf) - 5A{fip)xx 



(x, /ip) = xx{x - x)B{fil) 



(18) 

For = 3, in other words, for momentum transfers below the charm-quark production 
threshold, the functions A{fip) and B{fip) are defined by 



QB^, 



165? 





48 

Bl 




.as(/io). 
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"s(/io). 


."s(/io). 


48 
81 

+ 55f 


as(/^|) 

.tts(/io) 



101 

81 



101 
81 



(19) 



The DA of the octet rjs state contains only the quark component 0^g(x,/i|). This DA is 



identical to 0^^(a;, /i|), but with A{np) replaced by C(/i|), i.e. 



C(/^f) 



65? 



tts(/4) 

as(/io) 



50 
81 



(20) 



The explicit expressions for the functions A{^p) and -B(/i|) at momentum transfers above 
the charm quark threshold (or, for n/ = 4) can be found in the Appendix of Ref. 3|- For 
Uf = 4, the function C(/ip) should be modified to read 50/81— 2/3. If necessary, we shall 
distinguish between input parameters in Eqs. (|19|) and (|20|) by using the notations 51 (^71) 
and B^{r]s). 



III. BOREL RESUMMED 777 AND r]'-/ TRANSITION FORM FACTORS 

In Sec. ini we have outlined the key ingredients pertaining to both the standard HSA as 
well as the RC treatment of the transition FFs F^^^^Q^) and F^^i^Q"^). Let us now turn to 
a discussion of the main differences between these two approaches, starting with the choice 
of the scales fi^ and fip. It is evident that if a physical quantity can be factorized, like Eqs. 
(0) and ©, then the left-hand side (LHS) cannot depend on artificial intrinsic scales or on 
the particular renormalization and factorization schemes adopted. But at any finite order of 
QCD perturbation theory, truncation of the corresponding perturbative series will give rise 
to a dependence on the scales fip and fi"^, as well as on the factorization and renormalization 
scheme (for an in-depth discussion of these issues, we refer the reader to the second paper of 
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Ref. 0). Because higher-order corrections in perturbative QCD computations rule, 
large for both inclusive and exclusive processes, reliable theoretical predictions require an 
optimal scale-setting that minimizes higher-order corrections. Typically, the factorization 
scale enters the NLO contribution to the hard-scattering amplitude of meson transition 
or electromagnetic form factors in the form ~ ln((5^//x|), so that taking /i| equal to 
eliminates this term. But in order to analyze the sensitivity of our results to a chosen value 
of we shall perform all analytical computations for yUp 7^ Q^. 

The choice of the renormalization scale is somewhat subtler because this scale enters not 
only the NLO contribution, but also as the argument of the running strong coupling aj^\i^. 
To discuss this question, consider first the scale of the strong coupling. One effective method 
to solve this problem is the Brodsky-Lepage-Mackenzie (BLM) scale-setting procedure (26| . 
In this framework, a large part of the higher-order corrections — namely, those originating 
from the diagrams with quark "bubbles" insertions — can be absorbed into the scale of the 
QCD coupling constant. When utilizing this new scale one finds the NLO correction to 
be significantly reduced relative to its initial value. The generalization of the BLM proce- 
dure to all orders of perturbative QCD led to the invention of the RC method and the IR 
renormalon approach (for a review, see Ref. [l^]). In the case of inclusive processes, it was 
proven by explicit calculation that all-order resummation of diagrams with a chain of (quark) 
bubble insertions into the gluon line gives the same results as the calculation of one-loop 
Feynman diagrams for the quantity under consideration using the QCD running coupling 
at the vertices. Moreover, the IR renormalon approach in conjunction with the "ultraviolet 
dominance hypothesis" enables one to estimate higher-twist corrections to a wide range of 
inclusive processes. 

This approach was used for studying IR renormalon effects in exclusive processes 
as well. For instance, (— /9o«s/4vr)" corrections to the Brodsky-Lepage evolution kernel 
V\x,y\ajy(^'^^\ were computed in Ref. jH, ll^l and renormalon- chain contributions to the 
pseudoscalar meson DA and the 7r°7 transition FF were taken into account in ji^. Similar 



investigations along this line of thought were performed in Refs. |3Ct 131 . 

In addition to loop-integration ambiguities, exclusive processes may receive power- 
behaved contributions from the end point regions due to the integration in a process ampli- 
tude over the longitudinal momentum fractions of the involved partons. In fact, in order to 
reduce the NLO correction, for example, to the pion electromagnetic FF ^^^(Q^), the renor- 
malization scale /i^ should be set equal to the typical four-momentum, flowing through hard 
gluon lines in the partonic subprocess q(f -|- 7* — gg' |26|. Choosing the scale this way, 
inevitably leads to a dependence on the longitudinal momentum fractions carried by the 
hadron's constituents. In the case of Ft^[Q'^\ the NLO contribution to the hard-scattering 
amplitude T^{x,y,Q^) contains a logarithm of the form \n (xyQ"^ / fi"^) , with x and y being, 
respectively, the longitudinal momentum fractions of the quarks in the initial and the final 
pion. Hence, the natural choice to eliminate this term would be to set n'^ = xyQ'^. But 
due to the convolution of the hard-scattering with the soft components (cf. e.g., Eq. Q), 
integrations over x and y appear that give rise to power corrections when approaching the 
end point x — >• 0, 1; y ^ 0,1 regions. Renormalizing the process amplitude at a scale close 
to the large external momentum makes such contributions less pronounced but at the 
expense of large NLO logarithms. Therefore, if we are to optimize our theoretical calcula- 
tion, we have to minimize NLO contributions while keeping under control power corrections 
in the end point regions. 
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Specifically, for the meson-photon transition we have 



^i = Q'x, Jii = Q'x, (21) 

because in the corresponding two, leading-order, Feynman diagrams the absolute values of 
the square of the momenta flowing through virtual quark lines are determined exactly by 
these expressions. In the standard HSA one "freezes" the scale of the QCD coupling constant 
/^R (T^r); by replacing x by its mean value < x >= 1/2 and performs then the integrations 
in Eqs. (0) and © with as{Q'^/2). Let us stay within the HSA and concentrate on the NLO 
corrections to the quark component of Eq. (jSJ. Omitting an unimportant, in the present 
context, constant factor, we get 



2aJ^]tix,iJ,l)0<PUx,f^l), (22) 



2 



where the function t{x, /i|) is 



X 



2 xlnx 



In ^x — + (3 + 2 In x) a - 9 



X 



(23) 



and a = ln(Q^//Xp). In deriving Eq. (j221) we used the symmetry property of the quark 
component of the rji state DA, valid also for the function 0^g(x, /i|), 

0^1 (a;, /^f) = {x, l4) , 0»78 = (Pvs ^f) • (24) 

The generalization of our analysis to encompass the gluon component of the FF is straight- 
forward. 

Applying the RC method, the same quark component of the 7717 transition FF takes the 
form 



2a,iQ''x)t{x,fil)(g)(l)lix,fil). (25) 

M7I 



After a similar analysis for the gluon component of the form factor Q"^ F^^^{Q^) , using the 
RC method, we find 

Q'i^M,(Q')f °° ~ 2a,iQ'x)g{x,fxl) ® 0^,,(x,/x^) , (26) 

with the function g{x,fip) being given by the expression 

g(x,al)= _ + 6-= ln3: + 2a^^ (27) 

X \ xj X 

by making use of the antisymmetry of the gluon DA 0^^(a;, /i|) under the exchange x ^x 

0^,,(a;,/x|) = -0^,,(x,/x|). (28) 
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The gluon component in the standard HSA has the same form ()26|) with the argument of 
as being replaced by 

Summing up, we can write the transition FFs Q'^F^^{Q'^) and Q'^F^^{Q'^) in the context 
of the RC method as follows 



+ 



27T 



and 



T^ n(x) (g) (f)n^{x,fil) + —as{Q'^x)t{x,iJ,l) (l)n^{x, fi^) 

Ztx 



(29) 



(30) 



But the integrations over x in Eqs. and (jHUI) . when retaining the x dependence of the 
QCD coupling as{Q'^x) [as{Q'^x)]^ lead to divergent integrals because the running coupling 
as(Q^x) [Q;s(Q^;r)] suffers from an infrared singularity in the limit x — [x — > 1]. This 
means that in order to perform calculations with the running coupling, some procedure for 
its regularization in the end point a; ^ 0, 1 regions has to be adopted. 

As a first step in this direction, we express the running coupling as{Q'^x) in terms of 
asiQ'^), employing the renormalization-group equation, to find [32[, 



otsiQ'^x) 



1 + Inx/t 



as{Q^)f3i ln[l + lnx/t] 
27r/3o 1 + lnx/t 



(31) 

with I3q and (3i 



where as(Q^) is the one-loop QCD coupling, t = An / PoasiQ"^) = ln((5^/A^) 
being the one- and two-loop coefficients of the QCD beta function 

2 19 
/5o = ll-2%, A = 51-yn/, 

respectively. Equation (j3ip expresses as{Q'^x) in terms of as{Q^) to the ~ ^^(Q^) order 
accuracy. 

Inserting ()3H1 into the formulas for the transition FFs ()29|1 and ()3()|1 . we obtain integrals 
which are still divergent, but can be computed using existing methods. One of them, applied 
in for the calculation of the electromagnetic pion form factor, allows one to obtain 
the quantity under consideration as a perturbative series in asiQ"^) with factorially growing 
coefficients C„ ~ (n — 1)!. Similar series may be found also for the transition FFs Q'^Fl^^{Q'^) 

as(g' 



oo 

E 

n=l 



477 



(32) 



But a perturbative QCD series with factorially growing coefficients is a signal for the IR 
renormalon nature of the divergences in ()32j) . The convergence radius of the series fl32|) is 
zero and its resummation should be performed by employing the Borel integral technique. 
First, one has to find the Borel transform BlQ"^ F^^ ]{u) of the corresponding series 13: 



B 



E 

n=l 



U 



n-1 



'n-l)\ 



(33) 
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Because the coefficients of the series ()32p behave hke C„ ~ (n — 1)!, the Borel transform 
()33|) contains poles located on the positive u axis of the Borel plane. In other words, the 
divergence of the series (jH^ has been transformed into pole singularities of the function 
B[Q'^Fli^]{u). These poles are exactly the IR renormalon poles. 

Now in order to define the sum (j32|) . or to find the resummed expression for the form 
factors, one has to invert B[Q'^Fl^^]{u) to get 



P.V. 



duexp 



Attu 



B 



2 I7i 

M7 



(34) 



and remove the IR renormalon divergences by the principal value prescription. These inter- 
mediate steps can be bypassed by introducing the inverse Laplace transforms of the functions 
in (jnU, i.e., 

1 1 r°° 

duexp[-u{t + z)]u'"'^, Reu > (35) 



and 



{t + zY r(z/ 

ln[t + z] 

JtTW 



/ duexp[—u{t + z)]{l—'yE — \nu)u. 
Jo 



(36) 



where T{z) is the Gamma function, 7^ ~ 0.577216 is the Euler constant, and z = Inx [or 
z=lnx]. 

Then using and ^ we find [lllisj 



Ayr roc 

a,{Q^x) = — / due-'''R{u,t)x-'' 
Pq Jo 



A 



Here, the function R{u, t) is defined as 



2/3i 

R{u, t) = 1 — ^""(1 — — In t — In m) . 

Po 



(37) 



(38) 



Having used Eqs. (EZ|), and (|nZ|) in Eqs. (EH) and (jnO)), and performing the 

integrations over x, we obtain the FFs Q^F^^{Q'^) and Q"^ F^^i^Q"^) within the RC method; 
viz., 

Q'Fl,^{Q') = fLNi 6 + A{fil) + (1 + A{fil)) due--'R{u,t)Q,{u)- 



5A(/i|) 



rfue""*i?(u,t)Q2(M) 



=^B{fil) / due-'''R{u,t)G{u) 
Po Jo 



(39) 



and 



Q'F'iQ'') = f^Ns { 6 + C(/x^) + 



(1 + C(/i|)) due-'^'Riu, t)Qi{u)- 



due'''^R{u,t)Q2{u 



(40) 



The functions Qi{u), Q2{u), and G{u) have the expressions 
Qi{u) 
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^5(2, /3)i_, - ^5(1, /3)2-u + 2«^5(2, /3)i-u + 3(a - 3)5(1 -u,2), (41) 



Q2{u) = ^5(3, (3)2-u - ^5(2, /3)3-„ + 2a^i?(3, /?)2-„ + 3(a - 3)5(2 - 3) , (42) 



and 



G{u) = ^5(1, /3)4-. + 6^5(2, /5)3-„ - 4^5(1, /3)3-n + 2a^5(l, /5)4-„- 

5(2, /?)3_„ - 6-^5(3, /?)2-„ + 4-^5(2, /?)2_„ - 2a-^5(2, /?)3_. , (43) 



where the standard notation for the Beta function B{x,y) 

T{x)Tiy) 



B{x,y) 



has been employed. 

After some manipulations, the functions Qi{u), Q2{u) and G{u) can be recast into the 
more convenient forms 

^ , , 2 2 2a l + 2a ^ a-3 

= - /o .A3 - 71 3i + 77^— + 3 



(1-m)3 (2-m)3 (1-m)2 (2-m)2 (1-m)(2-m) 

^ / N 2 4 2 2a 

V2(m) = 77; - 7^ + 



(2-m)3 (3-m)3 (4-m)3 (2-m) 

l+4a l + 2a a-3 
+ 6- 



{3-uy (4-u)2 ' (2-m)(3-u)(4-m) ' 

_ ^ 2_ 2 _ 5-a 3-a 

^''^~(4-m)3 (3-m)3 + (2-n)2 (3-m)2+ (4 - m)2 ' ^ ^ 

One observes that the FFs given by (j39|) and (j4(jp contain a finite number of single, 
double, and triple poles located at the points = 1, 2, 3, 4. In other words, employing 
expression (j37p . we have transformed the end point a; singularities in Eqs. (j29|) and 
(l3(Jj) into (multiple) poles in the Borel plane u. These poles are the IR renormalon poles and 
consequently the integrals in Eqs. ()39p and ()4U|) are just the inverse Borel transformations 
fl34|) . in which the Borel transforms B[Q'^fI^^^]{u) of the NLO parts of the quark components 
and the gluon component of the scaled FFs are, up to constant factors, proportional to the 
functions 

R{u,t)Qi{u), R{u,t)Q2{u), R{u,t)G{u). 

As we have emphasized above, the IR renormalon divergences can be cured by employing 
the principal value prescription, which we adopt in this work to regularize divergent integrals. 
Therefore, the integrals over u in Eqs. (j39p and (j4Up are to be understood in the sense of 
the Cauchy principal value. Removing these divergences, Eqs. and become just 
the Borel resummed expressions [Q"^ Fl^^{Q'^)Y'^^ and [Q"^ F^^{Q'^)Y'^^ for these scaled FFs. 

It is known 0, that the IR renormalon pole located at the point Mq = ^ of the Borel 
plane corresponds to the power-suppressed correction ~ (1/Q^)", contained in the scaled 
FFs. To make the discussion of this question as transparent as possible, let us for the time 
being neglect the nonleading term ~ in (|HT|) and make the replacement R{u, t) — > 1 in 
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()37|) . Then, the integrals in the resummed FFs with multiple IR renormalon poles at uq = n 
can be easily expressed in terms of the integrals with a single IR renormalon pole at the 
same point uq = n (see, Eq. (jS^) below), so that our formulas (jH^ and ()40j) will contain the 
integrals 

/ e_aU ^ / ^^^g2^)^n-l^^ ^ i ^^^(g) ^ (45) 

Po Jo n — u Jo n 
where f2n{Q) are the moment integrals 

fpiQ) = 7^, f"^ dkkP-'a^iP) . (46) 
Jo 

The integrals /p(Q) were calculated before using the IR matching scheme: 



Q 



n=Q L 



as{Q ) 



2tcp 



[n\-T{n + l, p\niQ/fXj))], (47) 



where /i/ is the infrared matching scale and r(?T, + 1, z) is the incomplete Gamma function. 
In Eq. dHj) {/p(/U/)} are phenomenological parameters, representing the weighted average 
of as{k^) over the IR region < < fij, and act at the same time as infrared regulators of 
the right-hand side (RHS) of Eq. The first term on the RHS of Eq. ()47p is a power- 

suppressed contribution to fp{Q) and models the "soft" part of the moment integral. It 
cannot be calculated within PQCD, whereas the second term on the RHS of Eq. ()47|) is the 
perturbatively calculable part of the function fp{Q), representing its "hard" perturbative 
tail. In other words, the infrared matching scheme allows one to estimate power correc- 
tions to the moment integrals by explicitly pulling them out from the full expression, and 
introducing new nonperturbative parameters fp{^i)- The same moment integrals f2n{Q), 
computed in the framework of the RC method (LHS of Eq. (j45p ) . contain information on 
both their soft and the perturbative components. Indeed, numerical calculations demon- 
strate that the LHS of computed by employing the principal value prescription, and its 
RHS — found by means of (P7j) for p > 2 — practically coincide with each other. Therefore, 
we can state that the scaled and resummed FFs (jH^ and pUj) contain power corrections 
~ (l/Q^)", n = 1,2,3,4. Hence the usage in phenomenological applications of both the IR 
matching scheme and the RC method seems legitimate. In fact, both methods have been 
used to calculate the pion's electromagnetic FF |35||] and the vertex function Q'^F^ig*g*{Q'^ ,uj) 
[l3l | . But the RC method has an advantage over the IR matching scheme because it allows 
one to compute the functions fp{Q) without introducing the new nonperturbative param- 
eters Hi and fp{fii). Moreover, using this method, the parameters fp{f^i) themselves can 
be computed in good agreement with model calculations and available experimental data 

HI 113. 



The power corrections ~ (l/Q^)" are important in the region of moderate and change 
the behavior of the scaled and resummed FFs (jH^ and ()40|) as functions of significantly, 
both qualitatively and quantitatively. In the present work we have to deal only with a finite 
number of IR renormalon poles. Their number and location, in the case under consideration, 
depend on the DAs (fTH|) used in the calculations. The asymptotic DAs of the rji and 773 states 
lead to only two IR renormalon poles at Mq = 1; 2. Distribution amplitudes, which include 
higher-order Gegenbauer polynomials n > 2, may lead to a series of IR renormalon poles at 
Mo = 5, 6. . .. Note that at small momentum transfers, in each integral f2n{Q) associated 
with the pole uq = n, the soft part dominates. In the context of the IR matching scheme 
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the integral f2n{Q) at = /if even consists of just the soft contribution. Restricting 
our considerations to contributions arising only from the nearest to the origin m = IR 
renormalon poles (which are, of course, the dominant ones), entails two problems: first, it 
reduces the accuracy of the numerical results and second, one loses information on the DAs 
of the rji and 773 states. Therefore, for the self-consistent treatment of the FFs (j39|) and (|4Up 
we should take into account contributions coming from all IR renormalon poles. 

The principal value prescription, adopted here to regularize divergent integrals over u, 
generates itself power-suppressed ambiguities (uncertainties) 

where {$g(Q^)} are calculable functions entirely determined by the residues of the Borel 
transforms B[Q'^F^^^]{u) at the pole q = Uq and {A^^^} are arbitrary constants. Taking 
into account these ambiguities in Eqs. (jHIJjl and pUj) leads to a modification of the Borel 
resummed FFs, amounting to 

[Q'fH^ iQ')Y'' - [Q'F^'^ iQ')Y'' + (Q')]^"' • (48) 

In accordance with the "ultraviolet dominance hypothesis" , the uncertainty in Eq. pH|l will 
allow us to estimate power corrections to the scaled FFs stemming from sources other than 
the end point integrations. Indeed, by fitting the constants {Nq} to the experimental data, 
one can deduce some information concerning the magnitude of such corrections. 

It should be clear that regardless of the methods employed for the computation of the 
form factors, in the limit Q"^ ^ oo these must reach their asymptotic values. The im- 
portant problem to be clarified is then whether our resummed expressions [Q'^Fri^{Q'^)Y'^^, 
[Q'^Frfr^{Q'^)Y^^ lead in the limit — cxo to their corresponding well-known asymptotic 
forms. For the sake of simplicity, we restrict ourselves to the /i| = case. To answer the 
question posed above, we first explore the — > oo limits of the DAs 4>^-^^{x, Q'^), 0^^(x, Q"^), 
and (f)r]s{x,Q'^). Because in Eqs. (fT3jl and (fT^ the eigenvalues 7!^ < and their absolute 
values increase with n for all > 2, going to the asymptotic limit only the quark component 
of the rji state DA survives, evolving to its asymptotic form, whereas the DA of the gluon 
component M (x, Q^) in this limit vanishes, i.e.. 



yx, Q') Qxx, (x, Q') . 



The same arguments apply also to the DA of the r^g state, consisting only of the quark 
component 

0r,g(x,(5^) '^—^ 6XX . 

In our case this means that the following limits are fulfilled 

A{Q'), BiQ'), C{Q')'^^0. (49) 
Moreover, we have to take into account that in this limit the term ~ a 

2(Q2) in the expansion 

of asiQ'^x) in terms of as{Q^) has to be neglected The latter requirement is equivalent 

to the replacement 

/■oo roo 

/ e-'''R{u,t)du^ / e-'^'du. (50) 
Jo Jo 
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Then we obtain 



2 N 1 res 



6 iflN, + /8 AT, 



Po Jo 



(51) 



But this is not the final result because in the integral above t = ln{Q^/A^) and its oo 
limit has still to be computed. The integral 



due 



-ut 



1 



9 



9 



1 — u 2 — u 



(2-n)3 {2-uy 

can be expressed in terms of the logarithmic integral 

df 

U{x) = P.V. / P- 
Jo mt 

after performing the integration by parts of the first three terms to obtain 
°° e-"*c/M 1 In A ln2A/i(A") e'^'^du 



10 {n — m)^ 2n'^ 2n 2 A" ' Jo {n — uY 
where A = / K^. Employing the formula 13j 



n A" 



X" 



nlnx 



M 

E 

m=0 



M > 1 



and retaining in the expressions 



A^ 



In A 



/2(A") 



terms up to 0(1/ In A) order, we finally find 



i--«.(g^) 



l-f«s(Q') 



(52) 



(53) 



(54) 



(55) 



It is worth noting that numerical constants and terms ~ In A in Eq. appearing due to 
Eq. cancel out in the final result. 

The Eq. fl55p for M= rj and rj' supplies the asymptotic limit of the corresponding transi- 
tion EEs. These limits can be readily obtained within the standard HSA by means of the 
asymptotic DAs of the r/i and 773 states. Stated differently, by explicit computation we have 
proved that in the Q"^ ^ 00 limit the Borel resummed expressions (j39|) and (j4(jp lead to the 
well-known asymptotic forms of the F^^((5^) and F^'^(Q^) form factors. 
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FIG. 1: Predictions for the scaled form factors as functions of of the 777 (left panel) and 
(right panel) electromagnetic transition. For the solid curves the designation is B^ijix) = 0. The 
dashed lines correspond to iJfl^/i) = 10; for the dash-dotted curves we use -BfC^i) = l^- The data 



are taken from Ref. 11 1 



IV. EXTRACTING THE r/ AND ??' MESON DISTRIBUTION AMPLITUDES 

In this section we perform numerical computations of the Borel resummed 777 and ?7'7 
transition FFs^ in order to extract the 7] and rj' meson DAs from the CLEO data. We 
shall also compare our theoretical predictions with those obtained with the standard HSA 
[6|,|36|J, the aim being to reveal the role of power corrections at low-momentum-transfer in the 
exclusive process under consideration. In our calculations below we shall use the following 
values of A, /Xq and /i| 

A4 = 0.25 GeV, /i^ = 1 GeV^ /x^ = (56) 

and we shall employ both the one-angle scheme (jl(J|) and also the two-mixing-angles scheme 
(fT^ . Eqs. (fTI?|) and (^01) will be evaluated using the two- loop approximation for the QCD 
coupling as(<5^): 



/3oln(QVA2) 



2A lnln(gVA^ 
PI ln(gVA2) 



(57) 



The results shown in Figs. Q-E] — with the exception of Fig. [3 — are obtained within the 
ordinary octet-singlet mixing scheme. In Fig. ^the predictions for the 177 and ?7'7 FFs are 
presented for -82(^1) = B^ij]^) = 0.02 and various values of i?2('7i)- One appreciates that 
without the gluon contribution (i^K^^i) = 0) both FFs are slightly below the data points, 
especially Q'^Fj^/^{Q'^). But their deviations are not dramatic and to improve the agreement 
with the data, one has to include the contribution coming from the gluon component of the 



^ Notice that in this Section "FF" means the scaled form factors. 
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FIG. 2: The la area in the — B\ plane estimated within the RC method by comparing the 
CLEO data and the theoretical predictions for the resummed and scaled transition FFs r/7, rl^. 

rji'-f transition FF. The corresponding results are shown in Fig. ^ by broken hnes. These 
numerical calculations demonstrate that the gluonic contribution is important at relatively 
low values of the momentum-transfer Q^. From Fig. ^ it is clear that the gluonic contri- 
bution, arising from the r]i DA with i^K^^i) > 0, enhances the transition FFs Q'^Ff^^{Q'^) 
and Q'^Fjj/y{Q'^) in the region 1.5 GeV^ < < 12 GeV^ while reducing their magnitude at 




FIG. 3: The rjj (left) and rj'^ (right) scaled transition form factors as functions of Q^. The central 
solid curves are found using the values i?2(^i) = ^2(^8) = 0-05 and i^fC^i) = 17. The shaded 
areas demonstrate la regions for the transition FFs. 
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FIG. 4: The dependence of the 777 (left panel) and (right panel) scaled transition form factors 
on the values of the decay constants /i and fs- The octet-singlet mixing angle is 9p = —15.4°. 
The solid curves in both panels are calculated using fi/ fn = 1.17, fs/ fw = 1-26. The long-dashed 
curves correspond to fi/ fn = 1.17, fs/fw = 1-24 and the short-dashed one in the left panel are 
found by employing the values fi/f-K = 1.17, fs/fn = 1.30. The dash-dotted curves in both panels 
describe FFs obtained with fi/f-K = 1.20 and fs/ fn = 1.28. 

> 12 GeV^. This effect is sizeable for the r]''y transition FF relative to its counterpart for 
777, in particular, for larger values of -82(^71) and for smaller values of i^K^^i). The impact 
of the gluonic contribution on the 777, r]''-/ transition FFs is quite understandable, recalling 
that the physical r] and rj' states consist predominantly of the flavor SUf{3) octet rjs and 
singlet rji states, respectively, with the rjij transition FF comprising also a gluonic part. 
Therefore, the r]''-/ transition FF should be and is more sensitive to the gluonic part. 

These features of the 777 and r]''-/ transition FFs determine the la region for the allowed 
values of the Gegenbauer coefficients -BK^i) = B^ijis) and .Bf (^71), plotted in Fig. |21 In 
other words, the 777 and r]''-) transition FFs, computed in the context of the RC method by 
employing the model DAs with input parameters belonging to the shaded region in Fig. |21 
describe the CLEO data with a la accuracy. 

In Fig. 121 we plot the la areas for the 777 and 77'7 transition FFs. If we were to consider 
the 7/7 and 7/'7 transitions separately, these areas would be larger than those shown in Fig. 
01 For the 7/7 transition, the upper bound of the la region can be extended towards larger 
values of Q"^ F^^IQ"^) . For the rj''^ transition, the lower bound of the corresponding la region 
can be shifted towards lower values of Q'^Fjj'.ylQ'^). But their joint treatment leads to the 
picture drawn in Fig. El 

A major problem in extracting the values of theoretical parameters from the experimental 
data is their stability against uncertainties inherent in the theoretical expressions. In the case 
under consideration, expressions (jH), (jSHI), and (pn|) for the 777 and ?7'7 transition FFs depend 
on the factorization scale, on the QCD scale parameter A, the decay constants /i, /§, and 
the octet-singlet mixing angle 6p. As we have explained in Sec. 11111 the renormalization scale 
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FIG. 5: The dependence of the rj'j (left panel) and r]'^ (right panel) scaled transition form factors 
on the octet-singlet mixing angle 9p. The solid curves describe the "default" choice 9 = —15.4°. 
The dashed curves correspond to ^ = —16.4° and the dash-dotted ones to ^ = —14.4°. 

A'-R (Pr) context of the RC method is determined by the hard-scattering dynamics of 

the underlying partonic subprocess and is not a free parameter. Our analytical expressions 
for the transition FFs, calculated by keeping fip ^ Q"^, allow us to analyze the dependence of 
the extracted parameters Blirji), Blirjfi), and Blirji) on the factorization scale /Xp. We have 
performed the computation of the 777 and ?7'7 transition FFs using the values yUp = (5^/2 
and /i| = 2Q^ and found that our prediction for the la area (Fig. |2)) is absolutely stable 
against these variations. This means that the FFs determined by the input parameters from 
the la area in Fig. |2l by varying the factorization scale, remain within the corresponding la 
regions shown in Fig. 121 Stated differently, the variation of /Xp in the limits /ip G [1/2, 2]Q'^ 
does not change (shift, rotate) the la area in Fig. |21 On the contrary, the variation of the 
QCD scale parameter A modifies the la region in Fig.|21 The entailed modifications shift the 
region along both axes, retaining, however, its form stable. Thus, computations performed 
with A4 = 0.26 GeV result in the following shifts: along the Blirji) = Blij]^) axis: ~ 0.005, 
along the i^fl^i) axis: ~ 1. Hence, the modification of the la area is ~ 9% in the first 
and ~ 6% in the second direction, respectively, the percentages being given relative to the 
central values (see, Eq. (jU^ below). 

The response of the central curves (Fig. E)) on variations of the decay constants /i, /s, 
and such due to jthe octet-singlet mixing angle Op within corresponding phenomenologically 
allowed ranges j^, is demonstrated in Figs. 0] and |S1 It is remarkable that under these 
variations the central curves remain entirely within the associated la areas for the 177 and 
r/'7 transition FFs. It turns out that the 777 transition FF is more sensitive to the value of the 
decay constant /§ than the ?7'7 one. The results for the t]''-) transition FF obtained by varying 
the constant /§ G [1.24, 1.30]/,^ at fixed fi = 1.17/,^ practically coincide with each other. ^ 



^ This is the reason why in Fig.^b) the FF corresponding to the values /i = I.IT/tt, fs = l-SO/jr is not 
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FIG. 6: The rj'y (left panel) and r]'^ (right panel) scaled transition form factors as functions of Q^. 
All predictions have been obtained within the ordinary mixing scheme and using the initial input 
parameters and (|56|) . The broken lines denote the FFs with the uncertainties included via Eq. 
(|l8|) . and using the following values of {Nq}, g = 1,2,3,4 : 0.9 (dashed lines); -0.6 (dash-dotted 
lines) . 

On the contrary, the rj''^ transition FF demonstrates a rather strong dependence on the 
decay constant /i, whereas the 777 one is stable under such variations (cf. the short-dashed 
and dash-dotted curves, respectively, in Fig.Hl^a)). Our computations with 9p = — 15.4°±1° 
confirm the conclusion drawn in Ref. ^ that the FF for the 777 transition is more sensitive 
to 6p than the one for the ?7'7 transition. Summing up, we can state that the modification 
of the central curves in Figs. El and El due to the variations of the decay constants and the 
mixing angle discussed above, does not exceed the level of ±3% of their values. 

In Sec. mil we have emphasized that the ambiguities produced by the principal value 
prescription, inherent in the RC method, affect the predictions for the transition FFs in 
accordance with Eq. (jlH)) . The ambiguity [Q'^ F^^\Q'^)]'^'^^ depends on the rji and ris DAs 
and also on the constant {Ng}. In reality, however, for given DAs of the rji and rjs states, 
the available experimental information allows one to extract constraints on {Nq}. To effect 
the influence of such contributions, we show exemplarily in Fig. IHl predictions for the FFs 
with and without such ambiguities, utilizing the expansion coefficients i^K??!) = Bliris) = 
0.05, i?2 (^7i) = 17. We find that in order that the FFs remain within the corresponding la 
regions, the upper and lower bounds, respectively, for the constants {Nq\ are provided by 
the values {Nq = 0.9} and {Nq = —0.6}. Hence, the 777 transition FFs with the ambiguities 
included, corresponding to {Nq = —0.6} {{Nq = 0.9}) at < 4 GeV^, are larger (smaller) 
than the FFs without such corrections and are, in addition, smaller (larger) for > A GeV^. 
For the ?7'7 transition FF we observe, qualitatively, the same behavior, but with — 5 GeV^ 
as the transition momentum scale from the small to the large (and vice versa) regions. In 



displayed. 
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FIG. 7: The 777 (a) and rl^ (b) electromagnetic transition FFs vs. . The sohd hnes correspond 
to the ordinary octet-singlet mixing scheme with parameters B\{r]\) = i?2(%) =0.02 and i?2(??i) = 
18. The broken lines are obtained within the two-mixing angle scheme. The dashed lines describe 
the situation with the same parameters as the solid curves. The parameters for the dash-dotted 
curves are i?2(f?i) = -^Kr/s) = 0.15, 51(^1) = 

any case, the uncertainties do not exceed the level of ±11% of the corresponding FFs in the 
region ~ (1.73 — 2) GeV^ and reach a mere =f3% level in the region ~ (16 — 20) GeV^. 

Below, we present sample estimates for the eigenf unctions expansion coefficients of the 
rji and rjs DAs in the context of the ordinary mixing scheme: 

B^ir],) = , Bl{r]s) = , 5f (r^i) G [17, 21.5] , (58) 

Bl{r],) = 0.02 , Bl{7]s) = 0.02 , B'^{r],) e [16, 20.5] , (59) 

Bl{r],) = 0.05 , Bl{7]s) = 0.05 , 5f (r/i) G [15, 19] , (60) 

and 

51(771) = 0.1 , BliT^s) = 0.1 , Sf (^1) e [13.5, 17] . (61) 

The constraints (jKHjl - (jHT|) on the input parameter -BK??!) are extracted for fixed coefficients 
52(771) and B2{'r]s), and represent the la range for the values of Sfl^i) compatible with the 
CLEO data. Restrictions on the parameters i?2('7i) and B^iris) at fixed value of 52(^71) can 
also be derived. For example, for Blirji) = 16, we get 

S|(r7i) = 16, 5|(77i) = Bl{7]s) G [0.02, 0.11] . (62) 

Summarizing this point, the estimates for the Gegenbauer coefficients i?2(?7i), -BK^s) and 
Blirj^) in the DAs for the r/i, r/s states are 

Bl{r^i) = Bl{7]s) = 0.055 ± 0.065, 5|(r7i) = 18 ^ 4.5. (63) 
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FIG. 8: The r/7 scaled transition FF vs. . In the computations the ordinary octet-singlet mixing 
scheme is used. The upper (lower) bundle of curves is found within the standard HSA (RC method). 
The correspondence between the curves and the input parameters is as follows: for the solid curves 
B\{r]x) = B^{r]8) = 0, (r?i) = 0; for the dashed lines S|(r?i) = B^{r]s) = -0.05, (r/i) = 0; 
for the dash-dotted ones i?2(??i) = Blirjs) = 0.1, iJfC??!) = 0, and for the short-dashed curves 
S|(ryi) = Bl[^8) = 0.1, S|(r?i) = 15. 

Here some comments concerning the usual octet-singlet mixing scheme (fTUI) and the 
parameter set pi|) are in order. These parameters were extracted from the analysis of the 
CLEO data using the two-mixing-angles scheme, but staying within the context of the hard- 
scattering approach of perturbative QCD. Our computations demonstrate that adopting 
the RC method, the parameters given by (fTTjl satisfactorily describe these data, provided 
one uses the usual octet-singlet mixing scheme. Therefore, one can consider the parameters 
flll|) as a prediction of the RC method and the one-angle mixing scheme. This prediction 
differs from those obtained already within the one-angle mixing scheme, but employing the 
traditional theoretical methods (see, for example, Ref. [s^) 

However, our calculations do not exclude the usage of the two- mixing- angle scheme in 
conjunction with the RC method. But in such considerably larger contribution of 

the nonasymptotic terms to the DAs of the rji and rj^ states would be required. Carrying 
out such a computation via (fT^ . (fT^ . we obtained the results shown in Fig. [7| Inspection 
of Fig. El^a) reveals that the 777 transition FF found within this scheme lies significantly 
lower than the data. Therefore, to improve the results, a relatively large contribution of 
the first Gegenbauer polynomial to the DAs of the rji and rj^ states seems necessary. In 
Fig. [7|we display the FFs obtained using the parameters Blirji) = 0.15, -82(^78) = 0.15 and 
B^ivi) = 18- We consider the values -Bfl'?!) = ^tivs) = 0.15 as determining the lower 
bound for the admissible set of DAs in the context of the two-mixing-angles parametrization 
scheme. Hence, in the two- mixing- angles scheme, we obtain 

Blivi) = Blivs) = 0.15, i?f (r^i) G [16, 20] . (64) 

The 7] and rj' DAs were extracted from the CLEO data on the 177 and r]''-/ transition FFs 
having also recourse to the r/'-meson energy spectrum in the decay T(15') rj'X [36.]. In 
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FIG. 9: The ratio R{Q'^) for the 777 FF. The sohd hne corresponds to the input parameters -BfC^i) — 
-^K^) = -^2(^1) = 0- The dash-dotted curve describes the same ratio, but for i^fC^i) — ^K^) — 
0, i?2 (^1) = 14, while the dashed one corresponds to i?2(^i) = -^K^s) = 0.05, iJf (^1) = 10. 

both these papers the standard HSA was employed. In Ref . , estimates for the parameters 
BKfil), .Bfl/^o)' ^iif^o) were made within the two- mixing-angles scheme (fT^ . reading 

Bl{l GeV^) = -0.08 ± 0.04 , 5f (1 GeV^) = 9 ± 12 , B^{1 GeV^) = -0.04 ± 0.04 . (65) 

These coefficients are related to ours through the expressions 

B'M) = ^ , Bli^^l) = ^ , B|i^.l) = ^ . (66) 

Using our approach and the one-angle mixing scheme, the values of these parameters were 
determined to be 

Bl{l GeV^) = -0.214 ± 0.025 , 5f (1 GeV^) = 19.25 ± 2.25 , 5^(1 GeV^) = , 

^2^(1 GeV^) = -0.183 ±0.025, 5f(l GeV^) = 18.57 ±2.25, B^{1 GeV^) =0.02, 
^2^(1 GeV^) = -0.139 ± 0.022 , 5f (1 GeV^) = 17.8 ± 2 , 5^(1 GeV^) = 0.05 , (67) 
^2^(1 GeV^) = -0.07 ± 0.02 , 5f (1 GeV^) = 16.85 ± 1.75 , 5^(1 GeV^) = 0.1 , 

and 

^2^(1 GeV^) = -0.1128±0.045, 5f(l GeV^) = 16.208±0.144 , 5f(l GeV^) = 0.65±0.045. 

(68) 

In the case of the two-mixing-angles scheme, we find 

^2^(1 GeV^) = -0.050 ± 0.022 , 5f (1 GeV^) = 20.4 ± 2 , 5^(1 GeV^) = 0.15 . (69) 

One observes that within the two- mixing-angles scheme, the parameters B2, -Bf obey the 
constraints Bl{lGeV^) < and ^^IGeV^) > (cf. Eq. 
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On the other hand, the constraints for the parameters SK/^o) -^K/^o)' extracted in 
Ref. ||36j at the normahzation point /ig = 2 GeV^, read 

5|(2 GeV^) = 0.010 ± 0.068, 5f (2 GeV^) = 5.6 ± 3.4 . (70) 

Comparing now Eq. (f7n|) with the values given in Eq. ()64|) . and taking into account that in 
Ref. |3^ different values for the scheme parameters defined by Eq. were used, we come 
to the conclusion that in the context of the RC method and the two-mixing-angles scheme, 
the region -61(2 GeV^) < should be excluded as contradicting the CLEO data. 

The apparent discrepancy between the results of the present work and those of Ref. 
as regards the extracted values of the coefficients B^l GeV^) and -Bf(l GeV^), is 
related to the fact that the employed theoretical schemes are intrinsically different. Indeed, 
the transition FFs computed in the standard HSA overshoot the CLEO data — especially 
in the low-momentum transfer regime. In Fig. |S1 the 777 transition FF obtained in the 
standard HSA and the ordinary octet-singlet mixing scheme is depicted. One appreciates 
that the deviation from the data is considerable. The DAs corresponding to the parameters 
B2{r]i), Blirjs) > even increase this disagreement, whereas by adding the gluon component 
with i?2(?7i) > one can reduce it. Therefore to decrease the magnitude of the FFs, and 
achieve this way a better agreement with the data, the standard HSA would call for the 
two-mixing-angles scheme and for DAs mainly with i?2(?7i), -82(^3) < 0. The inclusion of 
power corrections changes this situation radically. In fact, at low- momentum transfer these 
corrections enhance the absolute value of the NLO correction to the FFs by more than a 
factor of 2.5 — 3 and, because the contribution of the NLO term to the FFs is negative, power 
corrections reduce the leading-order prediction for the FFs considerably, while at the highest 
values measured by the CLEO collaboration this influence becomes more moderate. As 
a result, the 777 and r]''-/ transition FFs computed using the input parameters from the la 
area in Fig. |21 within the RC method in conjunction with the one-angle mixing scheme are 
in agreement with the CLEO data. In order to quantify these statements, we show in Fig. 
inithe numerical results for the ratio 



Ires 



NLO r^t \ 

for some selected values of the expansion coefficients. 



V. CONCLUDING REMARKS 

In this work we have performed a computation of the 777 and ?7'7 transition FFs within 
the RC method. The latter has enabled us to estimate a class of power corrections to the 
FFs related to nonperturbative effects arising from the dependence of the strong coupling 
on the longitudinal momentum fractions of the partons inside the 77 and rj' mesons after the 
identification of the renormalization scale with a physical momentum depending on these 
fractions. This has been achieved by regularizing the infrared singularities ensuing from the 
end points x = 0, 1 by means of the principal value prescription within the IR renormalon 
approach. The effect of power-suppressed ambiguities to the considered form factors was 
addressed and their influence was found to be less important, though not negligible, with 
contributions varying in the range between 3% at high to 11% at low values. 

Contributions to the FFs from the valence quark as well as the two-gluon Fock-state 
of the 7] and 7]' meson DAs have been taken into account. We have obtained the Borel 
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resummed expressions [Q"^ Fu'y{Q'^)Y'^^ for the FFs and proved that in the asymptotic hmit 
Q"^ ^ CO they lead to the standard HSA predictions. 

We have demonstrated that the effect of the calculated power corrections on the rj'y and 
ri''y transition FFs is considerable. Indeed, at moderate values of the momentum-transfer 
< 5 GeV^ they turn out to enhance the absolute value of the 0{as) correction to 
the FFs more than 2.5 — 3 times. The ratio RwyiQ"^) of the corresponding contributions 
depends on the specific M7 transition under consideration and on the input parameters 
(Gegenbauer coefficients) of the rj and t]' meson DAs. These features of the power corrections 
have important consequences: the enhanced (negative) NLO correction significantly reduces 
the leading-order contribution to the FFs, so that the input parameters of the t] and t]' 
meson DAs, which correctly describe the CLEO data within the RC method, must obey 
the constraints presented in Fig. El by the shaded la area to fulfill Eq. It is worth 

emphasizing that our predictions for the f] and rj' meson DAs disagree with those extracted 
from the CLEO data in the context of the standard HSA. 

The DAs of the r] and t]' mesons obtained in this work can be useful in the investigation 
of other exclusive processes that involve t] and rj' mesons, especially at lower momentum- 
transfer values, where the standard HSA is most unreliable. 
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